Stability of the Einstein Static Universe in open cosmological models 
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The stability properties of the Einstein Static solution of General Relativity are altered when 
corrective terms arising from modification of the underlying gravitational theory appear in the cos- 
mological equations. In this paper the existence and stability of static solutions are considered in the 
framework of two recently proposed quantum gravity models. The previously known analysis of the 
Einstein Static solutions in the semiclassical regime of Loop Quantum Cosmology with modifications 
to the gravitational sector is extended to open cosmological models where a static neutrally stable 
solution is found. A similar analysis is also performed in the framework of Hoi'ava-Lifshitz gravity 
under detailed balance and projectability conditions. In the case of open cosmological models the 
two solutions found can be either unstable or neutrally stable according to the admitted values of 
the parameters. 
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I. INTRODUCTION 

The Einstein Static (ES) Universe is an exact solution 
of Einstein's equations describing a closed Friedmann- 
Robertson- Walker model sourced by a perfect fluid and 
a cosmological constant (see, for example, [Ij). This solu- 
tion is unstable to homogeneous perturbations as shown 
by Eddington f2i, furthermore it is always neutrally sta- 
ble against small inhomogeneous vector and tensor per- 
turbations and neutrally stable against adiabatic scalar 
density inhomogeneities with high enough sound speed 

In recent years there has been renewed interest in the 
ES Universe because of its relevance for the Emergent 
Universe scenario [Ij in which the ES solution plays a 
crucial role, being an initial state for a past-eternal in- 
flationary cosmological model. In the Emergent Universe 
scenario the horizon problem is solved before inflation be- 
gins, there is no singularity, no exotic physics is involved, 
and the quantum gravity regime can even be avoided. 
This model, relying on the choice of a particular initial 
state, suffers from a fine-tuning problem which is amelio- 
rated when modifications to the cosmological equations 
arise but then a mechanism is needed to trigger the ex- 
panding phase of the Universe (see [5]). 

The existence of ES solutions along with their stabil- 
ity properties has been widely investigated in the frame- 
work of General Relativity for several kinds of matter elds 
sources ([5] and references therein). ES solutions also ex- 
ist in several modified gravity models [7] ranging from 
the Randall-Sundrum braneworld scenario [5] to Gauss- 
Bonnet modified gravity and f{R) theories [9]. The issue 
of the existence and stability of ES solutions has also been 
considered in the semiclassical regime of Loop Quantum 



* E-mail: 'canonico@s a.infn .it| 
tE-mail: parisi@sa.inin.it 



Cosmology (LQC), in either the case of correction to the 
matter sector [TU] or the case of correction to the grav- 
itational sector [TT]. Recently the same issue has been 
also considered in the framework of Hoi^ava-Lifshitz (HL) 
gravity [12 and IR modified Hoi'ava gravity [T51IT^ . 

When dealing with higher order modified cosmologi- 
cal equations, the existence of many new ES solutions 
is possible, whose stability properties, depending on the 
details of the single theory or family of theories taken 
into account, are substantially modified with respect to 
the classical ES solution of General Relativity (GR). 

Often in such analysis the case of closed [k = I) cosmo- 
logical models is the only one considered, neglecting the 
intriguing possibility of static solutions in open (fc = — 1) 
cosmological models. Here we point out that, due to 
the aforementioned corrections to the cosmological equa- 
tions, open ES models may be found even in the case 
of a vanishing cosmological constant or when the perfect 
fluid has vanishing energy density. Examples are explic- 
itly provided for the case of LQC and HL gravity. 

This paper is structured as follows. In Sec. II, we con- 
sider static solutions in the framework of LQC, following 
and enlarging the analysis already performed in |llj . It 
is shown that, beside the ES solution of GR, a LQC solu- 
tion arises also in the case of open cosmological models; 
its stability is completely characterized. Following the 
same approach, in Sec. HI we consider static cosmolog- 
ical solutions in the context of HL gravity with detailed 
balance and projectability condition. Two solutions are 
found along with their stability properties. In Sec. IV, 
some conclusions are eventually drawn. 



II. LOOP QUANTUM COSMOLOGY 

In Loop Quantum Cosmology the quantization 
techniques borrowed by Loop Quantum Gravity, a 
background-independent nonperturbative quantum the- 
ory of gravity, are applied to symmetry reduced models 



(see [IS] and references therein). 

For the sake of simphcity, in this section we consider 
the modified Friedmann equations arising in the semi- 
classical regime of LQC [T^. We consider gravitational 
modifications only, neglecting the inverse volume correc- 
tion to the matter sector. The motivation is twofold: the 
analysis of this system allows a more transparent com- 
parison with the case of GR; moreover it allows us to 
follow the notations introduced in [TT] which will also be 
easily used in the analysis of the HL gravity presented in 
the next section. 

The model considered is sourced by a perfect fluid with 
linear equation of state p = wp plus a cosmological con- 
stant A. The classical energy conservation equation still 
holds, 
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while the loop quantum effects lead to a modification to 
the classical Friedmann equation, 



and to the Raychaudhuri equation 
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Notice that we are considering at once the fc = case 
and the fc = ±1 cases [TB]. Here k = 8ttG — Sir/Mp, and 
the critical LQC energy density is pc 
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A. Static solutions 

The system of Eqs.([T])-([3]) admits two static solutions, 
i.e. solutions characterized hy a = H = p = 0. The 
first solution corresponds to the standard ES Universe in 
GR; the second solution arises from the LQC corrective 
terms: 
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The conditions under which these static solutions exist 
are summarized in Table |lj they follow from a^ > and 
p > 0. The presence of the curvature index k is worth 
stressing, indeed the previous analysis [11] can be en- 
larged to enclose the fc = — 1 case where the two solutions 
still exist. 



B. Stability analysis 

The stability of the solutions Eqs.Q and (Is]) can be 
characterized using dynamical system theory and per- 
forming a linearized stability analysis. To this aim, we 



first have to rewrite the system of Eqs.Q-lls]) in the form 
of a genuine dynamical system. Indeed, in these equa- 
tions the three variables a, H and p appear but the actual 
dynamics is constrained on a two-dimensional surface de- 
scribed by the modified Friedmann equation. Thus, fol- 
lowing [11], we solve Eq.Q for 1/a^. Two solutions are 
found: 
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Substituting Eq.([6|) into Eq.Q, we find two branches 
for the time derivative of the Hubble parameter, thus 
the original system splits in a pair of two-dimensional 
nonlinear dynamical systems in the variables p and H: 

GR: p= -3Hp{l 
LQ:p = -3i/p(l 
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Each one of the systems (Is]) and (|9| admits a fixed point 
representing a static solution, that is. 



GR: H ^0 and po 
LQ : H = {) and po = 
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respectively. Substituting these values of po in Eq.([2J) 
one gets exactly the values of the constant scale factor in 
terms of the parameters as in Eqs.Q and (Is]). 

Finally, to characterize the stability of the solutions 
Eqs.Q and ^ we evaluate the eigenvalues of the Jaco- 
bian matrix for the two systems Eqs.Q and Q at the 
fixed points Eqs.pTj) and ([l2| ) respectively. 

For the system in Eq.(l8|), we recover the usual prop- 
erties of the ES solution in GR. The eigenvalues of the 
linearized system at the fixed point are 
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In the case of positive curvature index fc = 1, these are 
either real with opposite signs for A > and w > —1/3 
- thus the fixed point is unstable (of the saddle type) - 
or purely imaginary for A < and —1 < w < —1/3, 
so the fixed point is a center. In the case of negative 
spatial curvature index fc = — 1 , these are again real with 
opposite signs for A < and w < —1, so the fixed point 
is unstable (of the saddle type). In Fig. [I] an example of 
the latter case is depicted. 
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FIG. 1. Dynamical behavior of the system around the GR 
fixed point for the case fc = — 1,A<0, w<— 1 with A = 
-100, TO = -2, K = 25.13274123. 



For the system Eq. (|9| the eigenvalues at the fixed 
point are 



\lq = ±\/(kPc- A)(l + w). 



(14) 



In the case of positive curvature index fc = 1, the LQ 
fixed point is either unstable (of the saddle kind), when 
npc > A and — 1 < w < —1/3, or a center for the lin- 
earized system, i.e. a neutrally stable fixed point, when 
Kpc < A and w > —1/3. In the case of negative spa- 
tial curvature index /c = — 1, the eigenvalues are purely 
imaginary for Kpc > A and w < — 1, so we have a center 
for the linearized system again. In the latter case the 
fixed point is nonhyperbolic thus the linearization theo- 
rem does not apply. Nevertheless a numerical integration 
of the fully nonlinear system Eq.(|9| for initial conditions 
near the fixed point confirms the result of the linearized 
stability analysis (see Fig.l2|. It's worth stressing that in 
open LQC models a stable ES solution exists in the case 
of positive values of the cosmological constant as long as 

A < Kpc- 

The results of the linearized stability analysis are sum- 
marized in Table [H 



III. HORAVA-LIFSHITZ GRAVITY 

The Hoi^ava-Lifshitz gravity [17 is a power-counting 
renormalizable theory of (3-|-l)-dimensional quantum 
gravity. In the ultraviolet limit, the theory has a Lifshitz- 
like anisotropic scaling between space and time charac- 
terized by the dynamical critical exponent z = 3. In the 
IR limit the theory flows to the rclativistic value z — 1. 



FIG. 2. Dynamical behavior of the system around the LQ 
fixed point for the case fc = — 1, A < Kpc, w < — 1 with 
A = 10, lu = -2, K = 25.13274123. 
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TABLE I. Existence conditions and stability conditions for 
the static solutions in Eqs.Q and (|5|. 



The effective speed of light c, the effective Newton con- 
stant G and the effective cosmological constant A of the 
low-energy theory, emerge from the relevant deforma- 
tions of the deeply nonrelativistic z = 3 theory which 
dominates at short distances [17j : 
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The first of the equations in( 15 ) imposes a relation among 
the parameters c, Aw and A; thus, in order to have a real 
emergent speed of light c, for A > 1/3 the cosmological 
constant has to be negative Aw- However, after an an- 
alytic continuation of the parameters (see [TH]), a real 
speed of light for A > 1/3 implies a positive cosmological 
constant Aw- Thus, mimicking the notation introduced 
in p2], we introduce a two- valued parameter e — ±1, in 
order to examine both the aforementioned cases at once. 
The HL cosmology has been systematically studied us- 
ing dynamical systems theory in [20) . it has also been 
investigated in [21' using conservation laws of mechan- 



ics. Here we consider static solutions of the cosmological 
equations for the HL gravity when both the detailed bal- 
ance condition and projectability condition hold. 

First we recast the modified Friedmann equations of 
[18j in a form which allows an easy comparison with the 
formerly considered case of LQC [^ 

The modified Friedmann equation reads 
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and the modified Raychaudhuri equation reads 
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The conservation equation for the energy density of the 
perfect fluid still holds unchanged: 
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Besides the overall factor 



on the right hand side 
of Eqs.(|16|) and (|17[), the modifications to the cosmolog- 
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leal equations of GR consist of the higher order terms 
oc fc^/Aa"* which become dominant at short distance 
scales and do not afi^ect the classical cosmological equa- 
tions in the case of fiat models. 



A. Static solutions 

It can be readily found, imposing th e con ditions a — 
H = p = 0, that the system of Eqs.(18)-(17l admits the 
following two static solutions: 
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The conditions under which these static solutions exist 
are summarized in Tables Ull and ITlll 

The presence of the curvature index k and the param- 
eter e in Eqs.(19l and pO] ) is worth being stressed; in- 
deed the analysis presented in [T^] can be enlarged to en- 
close the fc = — 1 case where new interesting possibilities 
arise. For instance a physically meaningful ES solution 
is present even in the case of vanishing energy density of 
the perfect fluid, i.e. Eq.(19). 



B. Stability analysis 

The stability analysis can be easily performed reduc- 
ing the original system to an actual two-dimensional au- 
tonomous dynamical system by making use of the Fried- 
mann constraint. In this case the simplest and most 



^ According to the definitions given in Sec. II, c = 1 and k - 
Eq.(|16| and Eq.(|17|l have been written accordingly. 
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straightforward choice is to eliminate the dependence on 



p from the other equations, being Eq.(16l linear in p. 



This allows us to describe the dynamics with just one set 
of equations. Indeed, solving Eq.(16) for p. 
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and substituting into Eq.(17) one gets a flrst order non- 
linear differential equation. 
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which, together with the definition of the Hubble param- 
eter, 



d = aH 



(23) 



provides a genuine two-dimensional autonomous dynam- 
ical system in the variables a and H. The system admits 
two fixed points with energy densities as in Eqs.(19| and 



(20); thus, to characterize the stability of these solutions. 



we evaluate the eigenvalues of the Jacobian matrix for the 
system Eqs.(22) and (23) at the fixed points correspond- 
ing to Eqs.(19) and (20) respectively. 



The eigenvalues at the fixed point HLl read 
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For all the admitted values of the parameters this is a 
pair of purely imaginary eigenvalues thus the fixed point 
is a center for the linearized system. The point is nonhy- 
perbolic so the linearized analysis may fail to be predic- 
tive at nonlinear order, nevertheless a numerical integra- 
tion proves that this fixed point is actually a center (see 
Fig.|3|. 

The results of the stability analysis for the fixed point 
HLl are summarized in Table HTl 
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TABLE II. Existence conditions and stability conditions for 
the static solution HLl. 



The eigenvalues at the fixed point HL2 read 
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According to the admitted values of the parameters this 
is either a pair of purely imaginary eigenvalues, so the 
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FIG. 3. Dynamical behavior of the system around the HLl 
fixed point for the case k = —1 with e = 1, A > 1/3, A < 0, 
w > 1/3. 



fixed point is a center for the linearized system, or a 
pair of real eigenvalues with opposite signs, so the fixed 
point is unstable (of the saddle type). In particular, the 
solution is a center for — 1 < w < 1/3 and is a saddle for 
w<— loru;>l/3 (for an example of the latter case see 
Fig. 11. 
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FIG. 4. Dynamical behavior of the system around the HL2 
fixed point for the case k — —1 with e = 1, A > 0, A < 0, 
w > 1/3. 

The results of the stability analysis for the fixed point 
HL2 are summarized in Table ITTIl 
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TABLE III. Existence conditions and stability conditions for 
the static solution HL2. 



IV. CONCLUSIONS 

We have considered the existence of static solutions in 
the framework of two recently proposed quantum gravity 
models, namely, LQC and HL gravity. We have shown 
that the inclusion of a negative curvature index fc = — 1 
enlarges the ranges of existence of the solutions affecting 
their stability properties thus providing new interesting 
results. The solutions found display stability conditions 
rather different from those of the corresponding solutions 
in closed models and from the stability properties of the 
standard ES solution of GR. 

In the case of LQC gravitational modifications to the 
Friedmann equations, a negative curvature index allows a 
neutrally stable static solution with A < Kpc and w < ~1, 
in contrast to the GR case. In particular the LQC static 
solution exists and is stable in the case of positive values 
of the cosmological constant as long as A < Kpc- 

In the case of HL gravity two static solutions are found. 
The inclusion of the negative curvature index leads to 
a static solution (HLl) with negative cosmological con- 
stant and vanishing energy density which is neutrally sta- 
ble against homogeneous perturbations. Furthermore, a 
negative curvature index allows a static solution {HL2) 
which can be either a saddle, for w < —1 and w > 1/3, 
or a center for —l<w< 1/3. 

As already observed in the frameworks of different 
modified models [3 [101 |TT] , the regime of infinite cycles 
about the center fixed points must be eventually broken 
in order to enter the current expanding universe phase. 
To this aim a further mechanism is needed, whose anal- 
ysis is beyond the scope of this paper. 
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